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Abstract 

Double porosity models for the liquid filtration in an absolutely 
rigid body is derived from homogenization theory. The governing 
equations of the fluid dynamics on the microscopic level consist of 
the Stokes system for a slightly compressible viscous fluid, occupying 
a crack - pore space. In turn, this domain is a union of two inde- 
pendent systems of cracks (fissures) and pores. We suppose that the 
dimensionless size 5 of pores depends on the dimensionless size e of 
cracks: 5 = e r with r > 1. The rigorous justification is fulfilled for ho- 
mogenization procedure as the dimensionless size of the cracks tends 
to zero, while the solid body is geometrically periodic. As a result, for 
the long-time process we derive the usual Darcy equations of filtra- 
tion for the liquid in cracks, while the liquid in pores is blocked and 
unmoved. For the short-time processes the homogenized system con- 
sists of acoustic equations, describing a two-velocity continuum with 
three independent parameters: the liquid velocity in pores, the liquid 
velocity in cracks and the common pressure. The proofs are based on 
the method of reiterated homogenization, suggested by G. Allaire and 
M. Briane. 
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Introduction 

The paper concerns a liquid motion in an absolutely rigid porous bogy. By 
the present moment there are a lot of different mathematical models, de- 
scribing this physical process. They take into account a geometry of a space, 
occupied by the liquid (liquid domain), and physical properties of the liquid 
and the solid components. Among different models the simplest one is Darcy 
equations of filtration 

v = -kVq + F, V •« = (), (0.1) 
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for the macroscopic velocity v and the pressure q of the liquid, when the solid 
skeleton is supposed to be an absolutely rigid body and the liquid domain 
is a pore space. For more complicate geometry, when the liquid domain is 
a union of system of pores and cracks, there are different phenomenological 
models (see, for example, Ref. [3], Ref. [8], Ref. [IB] , Ref. [H]). Note, that 
pores differ from cracks by its characteristic size: if l p is a characteristic size 
of pores and l c is a characteristic size of cracks, then l p <C l c . The well-known 
phenomenological double-porosity model, suggested by G. I. Barenblatt , Iu. 
P. Zheltov and I. N. Kochina [3], describes two- velocity liquid continuum in 
absolutely rigid body, where macroscopic velocity v p and pressure q p in pores 
and macroscopic velocity v c and pressure q c in cracks satisfy two different 
Darcy laws 

v p = -k p Vq p + F, v c = -k c Vq c + F, (0.2) 
and two continuity equations 

V- v p = J, V- v c = -J. (0.3) 

The model is completed by postulate 

J = P{q c -q P ), (3 = const. 

Scientific and practical value of mathematical models describing such 
complicate processes, is obvious. But their physical reliability is also very im- 
portant. Namely, we say, that the given phenomenological model is physically 
correct if it is one of basic models of continuum mechanics (as, for example, 
Stokes equations describing a slow motion of a viscous liquid, or Lame's equa- 
tions describing a motion of an elastic solid body) or asymptotically closed 
to some physically correct phenomenological model on the microscopic level 
(that is a model, obtained by homogenization of some model on the micro- 
scopic level, depending on the small parameter). In their fundamental paper 
R. Burridge and J. Keller [S] have used the very natural scheme to justify a 
physical correctness of the well - known in contemporary acoustics and fil- 
tration phenomenological model of poroelasticity, suggested by M. Biot [I]. 
As a model of the porous medium on the microscopic level authors have con- 
sidered the mathematical model, consisting of Stokes equations describing 
liquid motion in pores and cracks and Lame's equations, describing motion 
of a solid skeleton. The differential equations in the solid skeleton and in 
the liquid domain are completed by boundary conditions on the common 
boundary "liquid domain - solid skeleton" , which express a continuity of dis- 
placements and normal tensions. The suggested microscopic model is a basic 
one, because it follows from basic laws of continuum mechanics ( see also E. 
Sanchez - Palencia [15J). After scaling there appears a natural small param- 
eter 5 which is the pore characteristic size l p divided by the characteristic 
size L of the entire porous body: 5 = l p /L. The small parameter enters both 
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into coefficients of the differential equations, and in the geometry of the do- 
main in consideration. The homogenization (that is a finding of all limiting 
regimes as 5 \ 0) of this model is a model, asymptotically closed to the ba- 
sic model and, by definition, is a physically correct. But even this approach 
is too difficult to be realized, and some additional simplifying assumptions 
are necessary. In terms of geometrical properties of the medium, it is most 
expedient to simplify the problem by postulating that the porous structure is 
periodic with the period 5. Under this assumption R. Burridge and J. Keller, 
using a method of two-scale asymptotic expansion, have formally proved a 
physical correctness of M. Biot's model. That is, the homogenized model 
derived by authors coincides with M. Biot's phenomenological model. 

For the same geometry of the pore space (let call such a geometry as a 
single porosity geometry and corresponding mathematical model as a single 
porosity model) in absolutely rigid solid skeleton, when a liquid motion is 
described by the Stokes system 

a^Av - Vq + F = 0, V • v = 0, 

for dimensionless microscopic velocity v and dimensionless microscopic pres- 
sure q, where 

&tj, — A i 2^ 2 , A*i = const > 0, 

r is the characteristic time of the process, pf is the mean density of the 
liquid, g is the value of acceleration due to gravity, \x is the fluid viscosity, 
and F(x,t) is the given dimensionless vector of distributed mass forces, the 
rigorous justification of the physical correctness of Darcy's law has been done 
by L. Tatrar (see Appendix in Ref. [15]). Later a physical correctness of M. 
Biot's models, under the same assumptions on the geometry of a pore space 
as in Ref. [B], has been rigorously proved in Ref. [10J - Ref. [14] . 

For more complicate geometry of the absolutely rigid body, when the 
liquid domain is a crack - pore space (let call such a geometry as a double 
porosity geometry and corresponding mathematical model as a double poros- 
ity model), that is when a solid skeleton is perforated by system of pores and 
system of cracks, some attempts to derive macroscopic models, asymptoti- 
cally closed to some phenomenological models on the microscopic level, have 
been made by T. Arbogast et al [2J, A. Bourgeat et al [5] and Z. Chen [7]. 
Because the last two papers repeat ideas of the first one, let us briefly discuss 
the main idea in Ref. [2J. As an initial model on the microscopic level the 
authors have considered a periodic structure, consisting of "solid" blocks of 
the size e surrounded by the fluid. The solid component is assumed to be 
already homogenized: there is no pore space and the motion of the fluid in 
blocks is governed by usual Darcy equations of filtration. The motion of the 
fluid in crack space (the space between "solid" blocks) is described by some 
artificial system, similar to Darcy equations of filtration. There is no any 
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physical base, but from mathematical point of view, such a choice of equa- 
tions of fluid dynamics in cracks has a very solid base: it is impossible to 
find reasonable boundary conditions on the common boundary "solid" block- 
crack space, if the fluid dynamics is described by the Stokes equations. But 
there are reasonable boundary conditions, if the liquid motion is described 
by Darcy equations of filtration. Therefore, the final homogenized models in 
Ref. [2], Ref. [5] and in Ref. [7] are physically incorrect. 

In the present publication we follow the scheme, suggested by R. Burridge 
and J. Keller [6], and look for double porosity models in absolutely rigid 
porous body, derived from the basic model on the microscopic level. In other 
words, we look for physically correct macroscopic models. 

We start with a liquid domain, composed by a periodic system of pores 
with dimensionless size 6 and a periodic system of cracks with dimensionless 
size e, where 5 = e r . The liquid motion on the microscopic level is described 
by the Stokes system (see Ref. [9]) 

dv _ da 

a T — = ail Av-Vq + F, + a q V ■ v = 0, (0.4) 

for dimensionless microscopic velocity v and pressure q of the liquid, where 



2/i L c 2 

n = — -■ " - - 
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and c is a speed of sound in fluid. 

The case r = 1 corresponds to already studied situation of a simple pore 
space, and the case r > 1 corresponds to a real double porosity geometry. 

We assume that all dimensionless parameters depend on the small pa- 
rameter e and the (finite or infinite) limits exist: 

lima T (e) = r , limaJe) = /i , lim a q (e) = c 2 f , lim — 7 = /ii, lim-^=/i 2 . 

In what follows we deal with cases, when 

/i = 0, r < oo, < Cj < oo. 

The aim of any homogenization procedure is to find all possible limiting 
regimes as e \ 0. Of course, these regimes for the model (10.41) depend on 
criteria r and //i, which characterize different type of physical processes. We 
may roughly divide all these processes into two groups: long-time processes 
(filtration) and short-time processes (acoustics). It is well-known, that the 
characteristic time of the liquid filtration is about month, while the charac- 
teristic size of the entire porous body is about thousand meters. Therefore, 
we may assume that for filtration r = 0. The rest of processes we call 
acoustics and all these situations characterized by criterion tq > 0. 
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If we consider the case of simple geometry (r = 1) in filtration, then 
the homogenization procedure has a sense only if /x 2 > (see Ref. [TU]). 
Moreover, if // 2 = oo (extremely viscous liquid), then the unique limiting 
regime is an immobility (the liquid is blocked in pores ). This fact is a 
simple consequence of the Friedrichs-Poincare inequality. The same situation 
is repeated for the case r > 1 of double porosity geometry. As before, the 
homogenization procedure has a sense if and only if \i\ > 0. But this criterion 
automatically implies the equality /i 2 = oo. Therefore, due to the same 
Friedrichs-Poincare inequality the limiting regime for the liquid in pores is a 
rest state (that is the velocity v p in pores is identically zero). Next we use 
the method of reiterated homogenization, suggested by G. Allaire and M. 
BrianepQ and for the velocity of the liquid in cracks v c and liquid pressure 
q, when the crack space is connected, we derive the usual Darcy equations 
of filtration: 

1 1 da 

v e = B<p -(F--Vq), -ji + V-v c = 0, v p = 0. (0.5) 

For disconnected crack space (isolated cracks) the unique limiting regime is 
a rest state: v c = 0. 

Note, that under condition /x 2 = oo the limiting regime for the liquid 
in pores in the single porosity model is also a rest state, as in the double - 
porosity model. Therefore, the main point here is not in the geometry 
of the solid skeleton, but only in physical assumptions. 

It is absolutely clear, that there is an overflow from pores to cracks and 
vice versa for the double porosity geometry, and there is a very slow flow in 
pores for the single porosity geometry in the real physical processes. There- 
fore, the initial postulates of a basic model (more precisely, the supposition 
that a solid skeleton is an absolutely rigid body) do not correspond to the 
reality. But, on the other hand, we can say that the obtained physically 
correct model describes the physical problem with a simulation error o(cJ ), 
where c s is one of the speeds of sound in the solid skeleton. 

Anyway, the double porosity models, suggested in Ref. [3], Ref. [2], Ref. [5], 
or in Ref. [7] have no common with a physically correct mathematical models. 
First of all, the correct models contain only one pressure, common both for 
pores and for cracks. Secondly, the liquid in pores is blocked and unmoved. 

For short-time processes (tq > 0) the situation is different. Under this 
restriction the homogenization procedure has a sense for any fix and /i 2 . In 
particular, for = and /x 2 = we show that the limiting regime is a two- 
velocity continuum with three independent characteristics: velocity in pores 
v p , velocity in cracks v c and common pressure q, which satisfy the acoustic 
equations, consisting of two momentum conservation laws in the form 

= K&I-Bf) ■ (--Vg + F), (0.6) 
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(m c I-Bi 2) ) • ( Vg + F), 



(0.7) 



and continuity equati 




c 



1 dq 



+ V ■ (v c + «p) = 0. 



(0.8) 



For sufficiently large cj (incompressible liquid) we may substitute (10.81) by 



which coincides with the sum of continuity equations in (10.31) . but equations 
( 10.61) and (10. 7p have completely different structure than equations (10. 2p . even 
for the case q p = q c . 

§1. Mathematical model on a microscopic level and basic 

a' priori estimates 

First of all we define the liquid domain Q £ , which is a subdomain of the 
unit cube Q. Let Q = Zf U Z s U 7 C , where Zf and Z s are open sets, the 
common boundary 7 C = dZf H dZ s is a Lipschitz continuous surface, and a 
periodic repetition of the domain Z s is a connected domain with a Lipschitz 
continuous boundary. The elementary cell Zf models a crack space Q e c : the 
domain Q e c is an intersection of the cube Q with a periodic repetition in M. 3 
of the elementary cell sZf. In the same way we define the pore space tip 
Q = YfUY s U7 P , 7 C is a Lipschitz continuous surface, a periodic repetition of 
the domain Y s is a connected domain with a Lipschitz continuous boundary, 
and Q 5 p is an intersection of Q\Q £ C with a periodic repetition in IR 3 of the 
elementary cell 5Yf. Finally, we put fl £ = fl p U fl e c . 

We may also characterize a liquid domain using indicator functions. In 
fact, let rj(x) be the indicator function of the domain Q in R 3 , that is r](x) = 1 
if x E f2 and i](x) = if x E IR 3 \fi. Let also Xp(y) De the 1-periodic 
extension of the indicator function of the domain Yf in Y and Xc( z ) be the 
f-periodic extension of the indicator function of the domain Zf in Z. Then 
Xc( x ) — v( x )Xc( x /£) stands for the indicator function of the domain Q £ c , 
Xp( x ) — 7 ?( a; )(l — Xc( x /£))Xp( x /$) stands for the indicator function of the 
domain fi^ and x £ ( x ) — X £ c ( x ) + Xp( x ) stands for the indicator function of 
the liquid domain Q £ . 

Next we define the generalized solution of the problem (10.41) as functions 
v £ and q £ such that 



V • (v c + u p ) = 0, 



v 



■ £ E L 2 ((0,T);W* (tt £ )) and 



o 



dq £ 
~dt 



E L 2 (n £ x (0,T)) 
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satisfying continuity equation in (10.41) in a usual sense almost everywhere in 
Q £ x (0, T) and integral identity 

J j (a T v £ ■ ^ - a^Vv £ : Vc^ + <f V • ip + F (p^jdxdt = (1.1) 

for any smooth vector-functions if, vanishing at dQ £ and t = T. The ho- 
mogeneous boundary condition v = on dfl £ is already included into the 

o 

functional space W\ (^ £ ) and homogeneous initial condition v(x,0) = is 
included into integral identity (II. ip . The known function F is supposed L 2 
integrable over domain Q x (0, T). In the same standard way, as in Ref. [TO] , 
one can show that for any e > there exists a unique generalized solution to 
the problem (10.41) and 

/ f (\v £ \ 2 +a^\X7v £ \ 2 + \q £ \ 2 + \V-v £ \ 2 )dxdt<c[ [ \F\ 2 dxdt, (1.2) 
Jo Jo? ^ ' Jo Jn s 

r (l£ l X £ p \v £ \ 2 dx + ^ [ x £ M\ 2 dx)dt<C [ T [ \F\ 2 dxdt, (1.3) 
Jo ^ o Jn £ Jn ' Jo Jn? 

where C is independent of e. 

To prove (11.21) and (11.31) we multiply the Stokes equation in (10.41) for v £ 

by v £ and integrate by parts over domain Q £ 

[ (a T \v £ \ 2 + — \q £ \ 2 )dx+ [ a^\Vv £ \ 2 dx= [ F-v £ dx. (1.4) 
Jn? ^ a q ' Jq s Jn? 

For r > (11.21) is a simple consequence of (11.41) . Let now r = 0. To estimate 
I s = L |f £ | 2 <ix we divide it by two integrals: 

i £ = i 5 p + i £ , i s p = [ \ v £ \ 2 dx, r c = [ \v £ \ 2 dx. 



1 d 
2di 



Let G p k \ where k = (ki,k 2 ,k 5 ) E Z 3 , be the intersection of with a set 
{x : x = e(y+ k), y e Y}. Then Q s p = U keZ 3G p k) and 

4 = E J ^)> = [j v£ \ 2dx - 

f /-Ate) 

In each integral I p we change variable by x = 5y, next apply the Friedrichs- 
Poincare inequality and finally return to original variables: 

\v £ \ 2 dx = 5 3 [ \v £ \ 2 dy < S 3 C {k) [ \V y v £ \ 2 dy = 5 2 C {h) [ | V x v £ \ 2 dx. 

Jy(fc) Jg^ 



7 



Here v £ (y,t) = v £ (x,t), Y^ k > C Y is an appropriate translation to origin of 
the set (l/6)G ( p k \ and is a constant in the Friedrichs-Poincare inequality 
for the domain Y^ k \ To estimate these constants uniformly with respect to 
5 (or e) let us clarify the structure of the domain Y^ k \ If the closure of 
Gp has no intersection with the boundary between pore and crack spaces, 
then Y^ = Yf and C^ k > coincides with a fixed constant C . Otherwise, 
y( fc ) is one of two domains, obtained after splitting Yf by some smooth 
surface, asymptotically closed to the plane as e \ 0. Due to supposition 
on the structure of the solid part Yf, constants uniformly bounded for 
all possible planes, splitting Yf. Therefore, supC^ fc ^ < C (for simplicity we 
denote all constants independent of e as C) and 

I S P <5 2 CJ2 f k \^ x v £ \ 2 dx<5 2 C ! \V x v £ \ 2 dx. (1.5) 

To explain ideas we consider the easiest geometry, when the liquid part Yf 
is "surrounded" by the solid part Y s . That is, for each facet S C dY of Y 
the liquid part S H dYf is completely surrounded by the solid part S H dY s . 
Next we extend he velocity v e from the liquid to the solid as zero. Then the 
constant in in the Friedrichs-Poincare inequality for Y^ depend only on the 
ratio a = Vf/V s between the volume Vf of the liquid part of Y^ k ' and the 
volume V s of the solid part of Y^: < Co. It is easy to see, that for 
chosen geometry of Yf and for any type of splitting of Y by planes, this ratio 
a is uniformly bounded. 

In the same way we show that 

I e c <e 2 C [ \V x v £ \ 2 dx. (1.6) 

Thus, 

s 2 f 2 r 

p<C(— + —)a a \V x v £ \ 2 dx. 

The rest of the proof is standard and follows from (jl.4p - (jl.6p and Holder 
and Gronwall's inequalities. 

§2. Homogenization procedure 

Now we ready to find limiting regimes for the model ( 11. 1ft . As the first 
step we extend the velocity and pressure from Q £ to Q as zero, keeping for 
simplicity the same notations. Next we use estimate (11.21) and conclude, that 
under condition 

r + /ii > (2.1) 
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(see Ref. [10]) the sequences {v 6 }, {q £ } and {V • v £ } weakly converge in 
L 2 (Q x (0, T)j (up to some subsequences) to v, q and V • v correspondingly 
and 

|| + V., = 0. (2.2) 

At the same time sequences {v £ } and {q e } three - scale converge (up to some 
subsequences) to 1-periodic in variables y and z functions V(x, t, y, z) and 
Q(x,t,y, z) correspondingly (see Ref. [1]) and 

Vdzdy = v, / / Qdzdy = q. 
'yjz JyJz 

Note, that three - scale convergence means the convergence of integrals 

v £ (x, t)-(p(x, t, — , %)dxdt — > / / / / V(x,t,y, z)-(p(x,t,y, z)dzdydxdt, 
to Jo JnJYJz 

for any smooth 1-periodic in y and z function <p(x, t, y, z). 

If we define a microscopic velocity in pores v £ p as v e = Xp v£ an d a mi- 
croscopic velocity in cracks v £ as v £ = x £ c vS then it is naturally to call a 
weak limit v p of the sequence {v £ p } as a macroscopic velocity in pores and 
a weak limit v c of the sequence {v £ c } as a macroscopic velocity in cracks. 
On the other hand, by definition of three - scale convergence, the sequences 
{v p } and {v £ c } three - scale converge to (1 — Xc( z ))Xp(v) V{ x i ^ Vi z ) an d 
Xc( z ) V{x, t, y, z) correspondingly and 




(1 - Xc)x P Vdzdy = v p , / / XcVdzdy = v c , Vp + v c = v. (2.3) 
z Jy Jz 

As the next step we state that Q = q(x, t)x(y, z)/m, where x — Xc(z) + (l — 
Xc{z))xc(y), rn = m c + (1 — m c )m p is the porosity of the crack - pore space, 
m c = J z Xc dz is the porosity of the crack space and m p = f Y x P dy is the 
porosity of the pore space. The proof of this fact repeats the corresponding 
proof in Ref. |10j . 

To find function V we must derive some microscopic periodic boundary- 
value problems in the domain W = {(y,z) : x{Vi z ) = !}• It is clear, 
that these problems somehow follow from the basic integral identity (11. ip as 
e \ 0. The condition suppy C Q £ for the test function ip in (11.11) means 
that cp = ip(x, t,x/e,x/5). This form of the test function creates difficulties 
in the limiting procedure. In following terms 

Jf(t)= I q £ V-ipdx= I q £ (V x -<p + -V z -(p + ):Vy-(p)dx, 
Jn Jn so 

J £ 2 (t) = I a^Vv £ : Vipdx = - / a^v £ ■ (A x¥ 5 + -V x ■ (V z <p) + 
Jn Jn v £ 
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gV x ■ (V y ip) + |rV z ■ (V y ip) + ^A,cp + ^A^v) dx. 

in the integral identity (11.11) appear unbounded items as e \ 0. 

The simple analysis of (11. 3p shows that for nontrivial homogenization 
procedure /ii < oo. Otherwise, v e p — > 0, v £ c — > as e \ 0. Next we separate 
different cases. 

2.1. Liquid nitration 

For filtration processes r = and H\ > 0. Therefore, /Z2 = oo and (11.31) 
implies i> p = 0. 

As a test function in (II. ip we choose the function cp = h(x,t)if)(x/s), 
where h(x,t) vanishes at dfl and at t — T and function t^(z) is finite in Zf 
(suppi/> C Zf) and solenoidal (V 2 - if) = 0). 

Due to relations 



lima T (e) = 0, lim — = 0, V y <p = 0, A y <p = 0, = V z -y> = 0, 

e\0 e\0 £ 

a three - scale limit in (11.11) yields integral identity 

III ((niV c ■ A z ip)h + (hF + —Vh) -ij>dz)dxdt = 0, 
Jo Jn Jz f x m J 

for the function 

V c = V(x,t,y,z)dy. 
Jy 

Reintegrating this identity with respect to variables (x, t) we arrive at 

0= / ( f i 1 V c -A z i/>+(F--Vq)-il>)dz= I ( f i 1 A z V c +F--Vq)-il>)dz, 
JZf v rn J J Zf m J 

and consequently 

liiA z V c {x, t, z) - V 2 n c (x, t, z) + F{x, t) - -Vq(x, t) = 0, ze Z f . (2.4) 

m 

The term V^n c appears due to condition V z ■ i/> = 0. 

Three-scale limit in the continuity equation in (10.41) in the form 



Jo Jn &t 



(2.5) 



with test function ip = £h(x,t)^o(x/e) yields 

V z ■ V c = 0, ze Z f . (2.6) 

Finally, the representation V c = Xc(z)V c and supposition [i\ > (for this 
case the function V Z (V C ) is L 2 -integrable in Zf) give us the boundary con- 
dition 

V c = 0, ze 7c . (2.7) 
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The problem (12TD . (I2^jl - flU]) is well-known (see Ref. p3] and Ref. [H)] 
and has a unique solution in the form 

V c = A«(z)-(F-ivg), 



or 



B«-(F--Vg), B« = / A«(z)(fe, (2.8) 
m J Zf 

which is exactly the Darcy law for the liquid in the cracks. Equation ( 12.81) 
together with continuity equation (12. 2p and relations (I2.3p . where v p = 0, 
constitute the homogenized double-porosity model: 

1 1 da 

v B = B® -(F-—Vq), ~^ + V-v c = 0, v p = 0. (2.9) 
m Cj at 

For disconnected crack space (we define a disconnected crack space as a struc- 
ture where the liquid part Zj is completely surrounded by the solid part Z s : 
dZ n 7 C = 0) the problem ([274]) . fl2T6|) - (J2IZD has a unique solution V c = 0. 
That is for disconnected crack space the liquid in cracks is also unmoved. 
2.2. Acoustics 

For acoustics tq > 0. As we have mentioned above, under this condition 
the homogenization procedure has a sense for any fii < oo. To explain ideas, 
we consider only the case fi\ = fi2 = 0. Under these restrictions the term 
If (i) in (11.11) goes to zero as e \ 0. It is clear, that we may pass to limit in 
If(t) only for test functions ip satisfying 

V z ■ if = V y ■ (p = 0. 

We have already used one of the possible ways to find microscopic equations 
(12.41) . To derive microscopic equations for pore space we need a test function 
if depending on both variable y and z. Such a choice is possible, if we put 

(p = h(x,t)\V z f c (-) x V y ^ p (-)J, supp ip c (z) C Z s , supply) C Yf. 

In fact, using the easy checking formula 

V • (a x b) = a ■ rot b — b ■ rot a 

we have 

V z • if = -h(Vyip p (y)) ■ (rot z (V z (p c (z))^J = 0, 
V y ■ <p = h(V z (p c (z)) ■ (xoty iy y <p P {y)) \ = 0. 



Therefore 



J!(t) = J q £ v x h- (y z cp c (j) x V^ p (|)) dx 
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Y JZ 



and the limit as e \ in (11. ip yields 

I I I I ( T o V ^+^'l V ^ h+hF )-( V ^c(z)xV y ^ p (y)yydzdxdt = 0. 

First, we reintegrate this integral identity with respect to variables (x,t) 

U (V z (f c (z) xV y (p p (y))dydz = 0, U = r — + —V x q - F, 

at m 

and after that, reintegrate the last identity with respect to variables (y,z): 

V y ■ (rotgU) = V z ■ (rot y U) = 0. 
This equation has a solution 

U = V y Q p (x, t, y, z) + Q c (x, t, y, z), 
1-periodic in variables (y, z), which leads to the desired microscopic equation 

dV 1 

T-o -jjr = - V y Q p - V z Q c V x q + F (2.10) 

at m 

in the domain W C Y <g> Z . Equation (I2.10p is completed with continuity 
equation (12. 6p and continuity equation 

V y -V = 0, (y,z)eW, (2.11) 

which is a result of three - scale limit in (10.41) for test functions in the form 



^ = 6h{x,t)^ Q {j)Mj)- 
Now we derive the microscopic equations in the crack space for functions 

V c (x,t, z) = J V(x,t,y,z)dy, U c (x,t,z) = J Q c (x,t,y, z)dy. 

To do that we just integrate (12.101) over Y and take into account the period- 
icity of the function Q p : 

T —^ = -V z Ii c V x q + F, zeZ f . (2.12) 

at m 

Equations (12.61) . (I2.12p are endowed with homogeneous initial and boundary 
conditions 

V c {x, 0, z) = 0, z € Z f , V c (x, t, z) ■ n(z) = 0,26 7c , (2.13) 
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where n(z) is the unit normal vector to 7 C at a point z G 7 C . The boundary 
condition in (I2.13P is a consequence of the representation V c = x c ( z )V c and 
continuity equation (12. 6p . 

The problem flZE]), (12121) . (12131) is well-known (see Ref. [15] and Ref. [H)]) 
and has a unique solution in the form 

V z U c = Af(z)-(--V x q + F), M^=[ kf\z)dz. (2.14) 
m J Zf 

Finally, to derive macroscopic equation (10.61) we just integrate (12.121) over 
domain Zf and take into account (12.141) . 

As a last step we must find some microscopic equations for two-scale 
microscopic velocity and pressure 

V p {x,t,y)= V(x,t,y,z)dz, U p (x,t,y) = Q p (x,t,y, z)dz 

J Z s J Z s 

in pores, when y G Y/. The main problem here is a definition of the term 
G(x,t,y) = Q c (x,t,y,z)n(z)dcr z , 

J 7c 

which appears after integration (I2.10p over domain Z s : 

r -kt = -Vj, n p - G + (1 - m c )( V, q + F), y eY f . (2.15) 

at m 

This function G stands for the interaction between pore and crack spaces. 
There is no any additional equation, which may define G, and to find at 
least the structure of homogenized equations for short-time processes in the 
model ( 10.41) we postulate the relation 

G = — -V x q + F), (3 = const > 0. (2.16) 
m 

Under restriction (I2.16P equation (I2.15p . continuity equation 

V tf .V p = 0, yeY f , (2.17) 
and initial and boundary conditions 

V p {x, 0, y) = 0, y G Y f , V p {x, t, y) ■ n(y) = 0, y G 7p , (2.18) 
uniquely define the homogenized equation 

ro^f = K&I-Bf ) " ( Vg + F), (2.19) 

at m 

where @ c = (3 + 1 — m c > (see representation (12.141) ) . 

Equations (I0.6p - (10. 8ft constitute the closed homogenized system, de- 
scribing acoustics in absolutely rigid body perforated by systems of pores 
and cracks, filled by slightly compressible liquid. 
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Conclusions 



We have shown how the new rigorous homogenization method, suggested 
by G. Allaire and M. Briane, can be used to clarify the structure of mathe- 
matical models for liquid filtration in natural reservoirs with very complicate 
geometry. Using the scheme, suggested by R. Burridge and J. Keller, we 
have shown that for double porosity geometry of the absolutely rigid 
solid skeleton 

1) the unique physically correct model for liquid filtration is a usual Darcy 
system for the liquid in cracks while the liquid in pores is blocked, 

and 

2) acoustics described by some physically correct model for two indepen- 
dent velocities in pores and cracks and for common pressure. 
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